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Abstract 

We propose a method for measuring the hadron charge asymmetry in e + e~ 
collisions which is based upon the fragmentation function formalism, and is 
largely independent of modelling of fragmentation effects. Furthermore, in this 
method, QCD radiative corrections can be accounted for in a systematic way. 
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1. Introduction 



Forward-backward asymmetries measured at e + e~ colliders near the Z resonance 
can be used to measure the electroweak parameter sin 2 #w/. The differential cross 
section for the production of a light fermion pair is given by 



^L- = a f (-(l + cos 2 6) + AL(s) cosO] , (1.1 



a cos a \8 

where 9 is the angle of the outgoing fermion / with respect to the direction of the 
incoming electron. The forward-backward asymmetry is defined as 

A FB = ^p^B (1-2) 

where a J (erf) is the cross section for producing the fermion in the forward (backward) 
hemisphere. At leading order on the Z peak (neglecting photon channel contributions) 
it is given by 

A' FB (Ml) = , (1.3) 

where 

v f = 4 - 2e f sin 2 6 W , a f = l( , (1.4) 

with ef and J 3 denoting respectively the charge (measured in units of the positron 
charge) and the third component of the weak isospin. 

Forward-backward asymmetries in the production of leptons are measured at LEP 
and SLD, and they constitute one very important measurement for electroweak tests 
(see for example refs. []T|). Due to the fact that Vf is very small for leptons, very 
high precision is required in this channel in order to get a useful measurement of 
sm 2 8]y. On the other hand, for quarks the vector coupling is not small. There, 
however, one has difficulty distinguishing between quarks of different flavours, or 
even between quarks and antiquarks. One approach is that of heavy flavour tagging 
in the hadronic final state (see refs. H§). Another possibility is to use instead the 
total charge asymmetry, averaged over all quark flavours, without flavour tagging. 
One defines 

' \U r_ ad 
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and one then needs to find a final state quantity which is related to the charge of the 
primary quark, at least on a statistical basis. 

Results from various approaches have been reported in refs. [4j-§|. The methods 
used there are all characterized by some "jet charge definition", which should be 
used to assign a charge to a jet. One then has to establish how often a jet with 
a given charge comes from a quark with the same charge. This is computed using 
Monte Carlo simulations. The drawback of this method is that it has to rely upon 
the hadronization models built into the simulations. A further drawback is that it 
is not at all clear how the strong corrections to the produced final state affect the 
measurement. One has no way to understand what part of the radiative correction is 
already incorporated in the simulation and what part is not. In fact, for a quantity 
defined in this way, it may be impossible to define the radiative corrections properly. 

In the present work, we propose a method for measurement of the charge asym- 
metry which overcomes some of the above drawbacks: it does not rely so heavily on 
hadronization models, and QCD radiative corrections are fully and reliably calcula- 
ble. Our method relies upon the reasonable assumption of valence dominance of the 
fragmentation functions at large values of the energy fraction x, which can be phrased 
in the following way: at very large x, only the hadrons containing a primary quark 
survive. 

In the following sections we give a detailed description of the method, illustrated 
by some Monte Carlo results obtained using the simulation program HERWIg! [0]. 

2. Fragmentation function formalism 

The single inclusive cross section for the production of a charged hadron is given 
by the formula 




1(1 + cos 2 9) 



da% 
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We define the total and asymmetric charged hadron fragmentation functions 




^We used the current version, HERWIG 5.7, with default parameter values. 
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F\ = /dcos0cos0 -— 7:-^-, ;r = — r- ■ (2.2) 

J \axa cos a ax a cos a J ax 

Following the notation of ref. || we have 

F h (x) = J2a ,f(s)(D>} + (x,s)+D<}-(x, S )) . (2.3) 
/ 

We use the annihilation scheme, in which the definition of the fragmentation function 
is such that the coefficient functions for F h have no radiative corrections if the scale 
in the fragmentation functions is set equal to the centre-of-mass energy squared, s. 
In the following we will drop the dependence on the scale, assuming that it is always 
set equal to s. 

In the annihilation scheme, the asymmetry coefficient function does receive higher 
order corrections. We have 

F h A {x) = Y.Ms){Df{x)-Df{x)) (2.4) 
/ 

where 

D){x) = D){x) - %^ f -(2 - z)D h f (x/z) . (2.5) 

This equation contains all that is needed to compute the radiative corrections of 
order a s . For the moment, we will consider only the leading order analysis, so we 
will assume that Dj = D'j. 

The electroweak coefficients (Xqj(s) and Af(s) are given by the formulae 



47TQ/ 2 

°o,/(s) = [ej + 2e f v e v f p 1 (s) + (v 2 e + a 2 e ) (vj + a'j) p 2 (s)} 



-If 2 

s 

„2 



A f( s ) = —— ea/[e/pi(s) + 2v e v f p 2 (s)] , (2.6) 



where 



1 sirrir 



4 sin 2 6 W cos 2 9 W (m| - s) 2 + m 2 z T 2 z 



P2(S) Usin^vKCOs 2 ^) (m| - s) 2 + m 2 z Y 2 z ' (2 ' ?) 
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On the Z peak, neglecting the photon s-channel contribution, we have 

/ \ 47TG! / 9 2\ / 2 2\ / \ 

a ,/(s) = — ^— (^ e +a e ){v f + a f )p 2 {s) 
^/( S ) = 2a e v e 2a,fVf P2\s) . [21 



3. Asymmetry from fragmentation functions 

Consider first the total and asymmetric fragmentation functions for pions. Ac- 
cording to our valence dominance hypothesis, for large enough x the tt + can only 
come from a u or a d quark. We then have 

F*{x) = a 0>u (s)DZ + (x)+a , d (s)Df(x) 

Fl{x) = A u (s) Dl + \x) - A d {s) Df (x) . (3.1) 



Using isospin symmetry and the on-peak expressions ( |2.8|) , we find 

Fl(x) ( 2v e a e \ 2a u v u - 2a d v d 



F*{x) W + vlJ (al + vl) + (a d + vj) 



0.0375 (3.2) 



assuming sin 2 9\y = 0.23, which we will use as our reference value. The asymmetry 
is small in the pion case because of the partial cancellation between the u and d 
contributions. 

To test the validity of valence dominance in a particular hadronization model, we 
studied the behaviour of the pion fragmentation functions using HERWIG. In fig. [I] 
we plot F\ (x) /F n (x) . The error bars correspond to 2.5 x 10 5 Monte Carlo events, i.e. 
about 5% of the current LEP data. The valence dominance value is also shown. We see 
that there is agreement for x > 0.5 at the 10 - 20% level. Corrections of this order are 
not unexpected, for example due to fast non-valence pions from resonance decays. It 
should be possible to estimate such moderate corrections from hadronization models 
without introducing a strong model dependence into the determination of sin 2 9w 
Recall that 2v e = 4 sin 2 9w — 1, so that a 10% precision in eq. (|3.2j ) implies about 1% 
accuracy in sin 2 9w- 

Let us now consider the production of protons. Again using valence dominance 
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Figure 1: Ratio of asymmetry and total fragmentation functions for charged pions, as 
a function of the energy fraction x. The dot-dashed line shows the valence dominance 
limit of eq. 



we get 



F p (x) 



a 0>u (s) [D?(x) + + a 04 (s) [D p d (x) + D*(x) + D p b (x)] 



Fl{x) = A u (s) [Dl{x) + DP(x)] + A d (s) [D>(x) + D>(x) + . (3.3) 



Observe that we have considered as valence dominant the production of protons from 
strange, charm and bottom quarks. In fact, the production of a strange, charm or 
bottom baryon will eventually lead through weak decays to a stable nucleon, and in 
no case (in the valence dominance regime) do we end up with an anti-baryon (baryon) 
from an initial quark (anti-quark). In the case of charm and bottom baryons, these 
decays will lead to a considerable degradation of the initial momentum. Since we 
intend to look at the large-x region, let us neglect for the moment the charm and 
bottom contributions. We then get 



x 



2v e a e 



Fp(x) 



2a u v u DP(x) + 2a d v d (D p d {x) + pr(x)) 

A + vl) K + vl)DZ(x) + (a d + vj)(D p d (x) + Dl{x)) 



(3.4) 



Unlike the case of light mesons, baryon production through resonances does not de- 
grade the baryon momentum substantially. Resonance production and decay tend to 
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randomize the charge of the final nucleon. We may therefore, as a first approxima- 
tion, assume equality of the up, down and strange fragmentation functions to produce 
protons at large x, to obtain 



Fp{x) 



0.1373 . 



(3.5) 



This value is much larger than in the pion case, since no cancellation takes place. In 
fig. H we plot F\(x) /F p (x). As one can see, the Monte Carlo results approach the 
value ( p.5|) very quickly. 
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Figure 2: As in fig. 0, but for protons and antiprotons. Dot-dashed: eq. ( |3~5l 



We can also obtain a larger asymmetry in the mesonic case by considering pi- 
ons and kaons together. We have in the valence dominance approximation (again 
neglecting the mesons coming from heavy flavour decays) 



2v e a e 



2a u v u D™(x) - 2a d v d (Df(x) + Df{x)) 



F*»{x) \a* + vl (a 2 u + v$D?{x) + (a d + vj)(Df(x) + Df(x)) 



(3.6) 



where m stands for both charged pions and charged kaons. In this case we would not 
expect equality of the three fragmentation functions involved. In fact, in the valence 
dominance limit the only charged meson into which a down quark can fragment is the 
7r~. By isospin symmetry, D% should equal . An up quark, on the other hand, 
can also go into a K + . We therefore expect D™ + > D™ . Our Monte Carlo studies 
suggest, however, that for x > 0.5 this difference is less than 20%. 
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In the fragmentation of strange quarks, one does not expect much leading pion 
production. There are however pions coming from neutral kaon decays, which en- 
hance the charged strange fragmentation function without contributing to the charge 



asymmetry, as assumed in formula (|3.6|) . These effects merit further study, both 
experimentally and theoretically. For now we will make the rough assumption that 
the three fragmentation functions are all equal to a first approximation. Under this 
assumption we obtain 



2v e a e 



2a u v u - Aa d v d 



F"(x) W + v 2 e J (al + vl) + 2(a d + vj) 



-0.078 



(3.7) 



i.e. about twice the value in the case of pions alone. Fig. |3| shows a comparison with 
Monte Carlo results. 
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Figure 3: As in fig. [I], but for charged pions plus kaons. Dot-dashed: eq. (|3.7|). 



For completeness, in Fig. [| we also show the total and asymmetric charge meson 
fragmentation functions F m (x) and F™(x). The values for F™ are divided by the 
limiting value of eq. (|3.7| ), to show the relation between the two functions at large 
x. The dashed curve shows the result of applying the higher order correction (|2.5| ) 
to the total fragmentation function. The effect is small compared with the probable 
magnitude of the corrections to valence dominance. 
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Figure 4: Monte Carlo results on the asymmetry and total fragmentation functions 
for charged mesons, (pions plus kaons), as functions of the energy fraction x, at 
E c .m. = 91.2 GeV. The asymmetry has been divided by its valence dominance limiting 
value of eq. ( ^.2| ). Dashed curve: result of applying the higher order correction ( |2.5| ) 
to the total fragmentation function. 

4. Conclusions 

The method we propose for studying hadronic charge asymmetries, and hence 
measuring sin 2 ^, is based directly on the fragmentation function formalism, for 
which the machinery of perturbative QCD is already well developed. In particular 
the relation between the total and asymmetric fragmentation functions for individual 
flavours is known, and measurements at different energies can be related, to next-to- 
leading order in a s . This makes the method attractive from the theoretical viewpoint. 

For phenomenological applications, one needs to exploit relationships between the 
fragmentation functions of quarks and their flavours. The type of valence dominance 
we have assumed seems less model dependent than other relationships between quark 
and jet properties that have been proposed for this purpose. Typical corrections are at 
the 10-20% level at high x and could be safely estimated from models. The predicted 
asymmetries for protons and light mesons at the Z resonance are substantial and of 
opposite sign, providing a valuable cross-check on the method. 

Other cross-checks may be obtained experimentally. One can, for example, require 
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one fast particle in one jet, and then study the fragmentation functions of the opposite 
jet. For example, if one requires a fast 7r + on one side, (which, according to valence 
dominance, can only come from a u or a d quark), this implies that on the opposite 
side we can measure the linear combination of fragmentation functions (v^ + al)D^ + 
(v% + a^)D^ directly, for any x value. This in turn can be used to test the valence 
dominance hypothesis itself, and to measure the K and p fragmentation functions for 
this combination. This measurement can in turn be used to correct for deviations from 
the valence dominance hypothesis. Many other combinations (for example requiring 
a fast proton or A s or a kaon on one side) can be studied, allowing for a wide variety 
of consistency checks and error estimates which are independent of hadronization 
models. 

We have concentrated here on light-quark hadrons, but of course the method 
can also be applied to heavy quark fragmentation. In particular the asymmetric and 
total D* fragmentation functions at large x would provide a rather clean measurement 
based on charmed quarks. 
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